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In this paper we consider the two-loop calculation of the disjoining pressure of a symmetric 
electrolytic soap film. We show that the disjoining pressure is finite when the loop expansion is 
£SJ . resummed using a cumulant expansion and requires no short distance cut-off. The loop expansion is 

resummed in terms of an expansion in g = Ib/Id where Id is the Debye length and Ib is the Bjerrum 
length. We show that there there is a non-analytic contribution of order gln(g). We also show that 
the two-loop correction is greater than the one-loop term at large film thicknesses suggesting a 
non-perturbative correction to the one-loop result in this limit. 
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I. INTRODUCTION 



The determination of the effective interaction between two surfaces of a film-like structure is essential for the 
understanding of the conformational stability of a whole range of systems in soft-condensed matter physics. Examples 
are found in interactions between membranes and colloid physics The effective interaction between two semi- 
infinite dielectrics separated by a vacuum was calculated by Lifshitz and the effect of an intervening dielectric 
was determined by Dzyaloshinski, Lifshitz and Pitaevskii Q and subsequently reformulated more simply in The 
effective interaction is due to van der Waals or dispersion forces and is given by the sum over the Matsubara modes 
of the problem; the non-zero frequency modes are a quantum effect and the contribution of zero modes corresponds 
I ' to a thermal Casimir type effect. In purely dielectric problems the contribution of each mode is described by a 
free Gaussian field theory and it has been shown how the input to these free field theories can be determined using 

■ dielectric data, see [§j and references therein. 

I/"") | In many situations, especially in biology, the middle slab of the system is filled with an electrolyte. For example, this 
is the case for a simple soap film connected to a bulk which is filled with an aqueous salt solution. It has been argued 
that the presence of an electrolyte will not effect the non-zero frequency contribution to the interaction basically 
because the response time of the ions is too large to couple them to the these non-zero frequency modes. However 
the zero-frequency, or static, component of the interaction does couple to the ionic distribution and consequently 
its contribution is strongly affected by the presence of electrolyte. The field theory describing the zero frequency 

■ fluctuations of the electric field in the presence of ions is no longer a free field theory; the interactions are made up of two 
components: the basic thermal fluctuations of the electrostatic field and the induced Coulombic interactions between 

I i the ions. In the presence of a neutral surface and a symmetric monovalent electrolyte the mean-field electrostatic 
potential is zero and the one-loop contribution about this trivial mean-field solution is equivalent to Debye-Hiickel 
theory. The one-loop result is described and computed in [j| and has been rederived using a variety of different 
methods, for instance see 0,0 ■ In the absence of electrolyte the zero frequency contribution to the effective interaction, 
commonly known as the disjoining pressure Pd, behaves as Pd ~ — A/l 3 for a film of thickness I. In the presence of 
electrolyte, at large film separations relative to the Debye length Id, the disjoining pressure is screened at one loop 
and has the form Pd « — Aexp(— 2l/ln)/l. In both of these cases the constant A is positive and hence the effective 
force between the surfaces is attractive at one-loop order. 

In this paper we calculate Pd to two-loop order. The naive two-loop calculation is in fact divergent when the 
dielectric constant of the film, or slab, is greater than that of the exterior but we show here that the result can be 
resummed using a cumulant expansion similar to that used in |10| for the bulk to give a finite result which we argue 
is to be expected on physical grounds. The expansion we use can be controlled systematically, and the perturbative 
parameter is g — Id/Ib where Ib is the Bjerrum length. While the one- loop or first term in the expansion in g is 
regular, the second order term in g is of the form g 2 \n(g), this singular behavior is found in the Onsager-Samaras 
limiting law for the excess surface tension of electrolyte solutions and has its origin in the image forces which repel the 
ions from the air/water interface. We also analyse the behavior of the 0(g 2 ) term for large film thickness I and find 
that this term has the behavior — A' \n(l/lr)) exp(— 2l/lr>)/l, which means that it becomes larger than the one-loop 
result for very thick films. We conjecture that this result suggests that it is the first term in a series which sums to 
give a contribution at large I to the disjoining pressure with behaviour Pd ~ — A" exp(— 2l/lr))/l a , where the exponent 
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a has the form a = 1 + a±g + dig 1 ■ 



II. TWO-LOOP CALCULATION 



The model of the soap film we shall consider is the following. The film body is a planar slab of thickness I surrounded 
on both sides by an external dielectric medium such as air. The interior of the film is filled with electrolyte and 
connected to a bulk reservoir of electrolyte. The electrolyte consists of a passive solvent of uniform dielectric constant 
containing monovalent anions and cations which are treated as point charges of charges — e and e respectively. The 
electrolyte is symmetric in the sense that the anions and cations are identical in every respect other than that they 
have opposite electric charges. 

In what follows we shall use the standard field theoretic formulation for symmetric monovalent electrolyte systems 
. The grand partition function for a film of electrolyte solution, surrounded on both sides by an external medium 
such as air, of thickness I is given by 



E F (l)= / <#]exp(S F [<M) , (1) 



where the action S F of the film is given by 



S F [(j>,t\ = ~ I dx (V0) z + 2/i / dx cos (e0<f>) - ^ / dx (V0) 2 . (2) 
2 J Fl Jfi 2 J Fe 

The region inside the film is denotes Fj and, if z is the coordinate perpendicular to the film surface, then Fj is the 
region z € [0,1]- Inside Fj the dielectric constant is e and is taken to be that of the solvent, in most cases water, and 
the fugacities for the cations and anions in this symmetric electrolyte are the same and are denoted by fx. The region 
outside the film is denoted by F F ; in F F there is no electrolyte and the dielectric constant is denoted by e . In this 
paper we consider the physically relevant case where e > eo- 

For a bulk solution of thickness I the grand partition function is given by 

3 B (/)= / d[cP]exp(S B [(t>,l}) , (3) 



where Sb is the bulk action 
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S B = ~ I dx (V0) 2 + 2(j, I dx cos (ejicj)) . (4) 



B JB 



Here there is no contact with an exterior region and one can close the system by using periodic boundary conditions 
in the z direction. The value of the fugacity is determined by the bulk density p and is given by 

/"* = 7 / P /ia\\ = Z P ■ ( 5 ) 

(cos (e/30)) B 

Here, Z^ 1 = (cos (e/30))s is a renormalization factor accounting for ion self-interactions, and the expectation is taken 
at an arbitrary position in an infinite bulk system, i.e. where / — > oo. In experiments on planar systems or films the 
physical quantity which is measured is the disjoining pressure which is due to the effective interaction induced by the 
presence of the film surfaces. The disjoining pressure is the excess film pressure Pp over the bulk pressure Pb and 
can be measured exactly using a pressure cell |l2j • 

P d {l)=P F (!)-P B ■ (6) 
The film pressure depends explicitly on the film thickness and is given by 

P F (l) = ~ln(E F (l)) . (7) 

The bulk pressure is defined in the thermodynamic limit and is given by 



3 



In order to develop a systematic expansion for the disjoining pressure we pass to the rescaled model in terms of the 
rcscalcd field 



ep 

and we measure length in terms of the Debye length Id 

x — * x l D , 



(9) 



(10) 



where m = y/2pe 2 f3/e is the Debye mass and Id = l/m. The dimensionless coupling g is given by g = Ib /Id where 
Ib = e 2 /3/Aire is the Bjerrum length. In units of the Debye length we denote the thickness of the film by L = I/Id- 
The disjoining pressure Pd is obtained from the disjoining pressure of the rescaled model Vd using 



P d = m 6 V d = 87rpgVd ■ 
The effective action for the film in terms of the rescaled fields and lengths is 

,2 , Z(g) 



s F = -I J pi rfx (w) 2 + m J pi rfx cos ( ^) - £ J fe & m f , 

where n = eo/e and the renormalization factor Z(g) simply becomes 

z{9) = (cos ym~ B ' 

where again the above expectation is for an infinite bulk system. The action Sf can be decomposed as 



S F = 



z(g) 

Ang 



V F +S 



(o) 



AS F 



(11) 



(12) 



(13) 



(14) 



where Vf is the volume of the film and the first term is the ideal contribution. The term S^p is the action for a free 
or Gaussian field theory and is given by 



S 



(0) 



1 

"8^ 



Fi 



(V^) 2 + 4> 2 



n 

8^ 



Fe 



The interacting part of Sp is expressed as a perturbation 



AS* 



1 

Ang 



f/x 



Fj 



Z{g) (cos (V?<^) - 1) + 



and the action Sb for the equivalent bulk system is given by 



<2x (V0) z 



zjg) 

Ang 



dx cos (y/~g~4>) ■ 



Clearly Sb is invariant under <f> — > —(f) and hence from Eq. (|13[) one must have the expansion 

Z(g) = 1 + z l9 + z 2 g 2 + z 3 g 3 + ... . 
Using the same decomposition for the bulk action as for the film action we obtain 



s B = ^E1 Vb — 
Aug 8tt 



rfx 



l 



rfx 



Z(g) (cos (y/g<t>) - 1) + 
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where Vb is the bulk volume. To order g the action Sb may thus be written as 



1 + zxg + z 2 g 2 

<->b — : ' b 

Ang 



-1/ d. 

8tt 



(V0) 2 + <f> 2 



Aw 



b 4 z^ 2 



4! 



0(9 2 ) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 
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Hence in order to calculate ln(Ss) to order g one needs to evaluate Z(g) to order g 2 . We define 

z n (g) = 



(21) 



where the notation (0) n signifies the expectation value of O evaluated to n— th order in the cumulant expansion. 
Then Z n is an approximation to Z correct up to and including 0(g n ). Then we can write Sb correct to order g as 



Z 2 (g) 1 f 

S B = —. V B - — dx 

Airg 8tt J b 



(V0) 2 + 2 



1 

Ang 



rfx 



Z 1 (g) (cos (V50) - 1) + *f- 



(22) 



Clearly, calculating with the above action gives the value of ln(Ss) correct to order g which includes the two- loop 
term in the loop expansion coming from the term gcj) 4 /4! in the interaction term. We might go further and expand the 
interaction term and keep only terms of order g. However, when we consider the computation of the partition function 
Sj? for the film we shall see that the naive expansion of the interaction term in this manner is illegal because it gives 
rise to a spurious divergence due to the image charge singularities. In the expansion for the disjoining pressure we 
must first identify the correct Boltzmann factor associated with the image charge potential, which is not expansible in 
g, and only then may we expand the remaining contributions to order g to obtain the two-loop divergence-free result 
for the disjoining pressure. We will find a term that behaves like gln<? which indicates that some contributions are 
indeed not expansible in g. Therefore, at two- loop order we use the form for Sb given in Eq. (|22|l and a similar form 
form for Sf which we discuss below. Then to O(g) we have 



where 



with 



and 



ln(S B ) = ^^V B + In(S fll0 ) + (AS) B ,o + 0(g 2 ) 
4ng 



2s,o = j d[4>] exp (S B .o) , 



W) 2 + 4? 



8tt Jb 



AS= / dx 

Aug 



Zi(g) (cos (V50) - 1) + 



(23) 



(24) 



(25) 



(26) 



The last term is just the first term in the cumulant expansion about the free field theory with action Sb and is given 
by 



(A5)b, 



/ d[(j>] exp (S B ,o) AS 



J d[<f>] exp (S b ,q) 
Inside the bulk we define 

<0(x)0(x)) B ,o = G B (0) , 
by translational invariance within the bulk. This now yields 

1 



Zi(g) 



(V^0))b,o 



exp 



(|G B (0) 



and 



V B 



ln(5 B ) = -t^Vb + ln(H B , Q ) + 

4vrg 47r<7 



+ |G B (0)-exp(|G B (0))] +0( 5 2 ). 



(27) 



(28) 



(29) 



(30) 



One now carries out exactly the same calculation in the film where within a film of thickness L we define 



(0(x)0(x))jr )O = G F {z,L) , 



(31) 
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where the field correlator at coinciding points now depends on the distance z from the surface of the film and L, its 
thickness; (0(x)</>(y))^ o is the propagator for the action iSp^. We then obtain to O(g) 

ln(S F ) = ?M VF + i n (E K0 ) + ^-J^d X [exp (-l(G F (z,L) - G B (0))) + 9 -G F {z,L) - exp (| G B (0))] + 0( 5 2 ) . 

In the limits L — > oo, z — > oo, z/L — > 0, we recover the expression for ln(Ss) in Eq. f|30|> - We note that the exponent 
of the first exponential under the integral sign in Eq. 132(1 contains the repulsive image charge potential which is 
singular as z — > 0. Evidently, it is illegal to expand this exponential as it stands and it is clear that it contains the 
correct Boltzmann factor for the potential due to the image charges. In terms of the rescaled system the disjoining 
pressure is then given by 



^ d = f L B (L) + f L B l{ L) 



where 



B (L) = - (m(£ F ,o(L))-£ r lirn 



ln(H B , (L')) 



(33) 



(34) 



is the one-loop or 0(1) contribution in the expansion in g, and A is the area of the film. The two-loop term, or term 
to 0(g), is given by 



Bi(L) 



Ang 



dz 



c^(- 9 -G R (z,L))-l+ 9 -G R (z,L) 



(35) 



with 



G R (z,L) = G F (z,L)-G B (0) 



(36) 



We note that there is no need to calculate ^2(5), the approximation to Z(g) accurate to 0(g 2 ), in order to compute 
the disjoining pressure, since this term is the same in both the film and the bulk and hence cancels identically. We 
find 



/oo 
dp 



2Y 2 (P) cxp(-2PL) + T(P) exp(-2P(L - z)) + T{P) exp(-2Pz) 



1 -T 2 (F) exp(-2PL) 



(37) 



where 



r(P) = p - n ^El . (38) 

V ; P + nVP 2 - 1 

We can compare the result in Eq. ((35(1 with the result of Netz |£| where he computes the two-loop contribution to the 
free energy of the film. When normalized with respect to the bulk free energy, the result of Netz agrees with Eq. I(35|) 
expanded to and truncated at O(g). In the case where A = the expression of Netz is indeed the disjoining pressure 
correct to two-loops. However, in the case where A = (1 — n)/(l + n) > 0, i.e. when there is dielectric discontinuity 
in the system, the expression truncated at O(g) is divergent. As already remarked, this is due to the presence of 
image charges which repel charge away from the surface into the film when A > as is the case of air/water films, for 
example. Netz obtains a finite result using this truncated expansion by introducing an artificial Ultra- Violet cut-off 
in the P integration defining Gr(z, L) in Eq. ((37(1 . The UV cut-off is associated with a microscopic length scale such 
as the ionic radius of the solute. Expanding the exponential in the expression for B\ amounts to an illegal expansion 
of a Boltzmann weight which keeps charge away from the surface. This term needs to be treated with care as in the 
Onsager- Samaras [l3| calculation of the excess surface tension of ionic solution where a similar term also arises. 
We write 

(iV d = f3V { p + /3V { d 1] , (39) 

where is the 0(g n ) contribution to Vd- The one-loop, O(l), contribution to the disjoining pressure has been 
worked out by many authors and is given by 

Pl d - 2W1 l-r2(P)exp(-2PL) • { U} 
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From above we find the two-loop contribution to be given by 



PPf = T~ ( ex P [~Gr{L/2,L)\ + 2-G R (L/2,L) - 1) + f- / dz^-G R (z,L) 1 — exp [~^G R (z, L) 



1 



Aug 



■in 



L/2 



d 



dL 



We emphasize that this expression is finite for all L > and requires no UV regularization. 
We may attempt to expand Eq. (|41|l to O(g) to get 



L/2 g 

dzG R (z,L)—G R (z,L) + 0(g 2 ) 



(41) 



(42) 



thus giving what naively looks like the two-loop contribution. 

The function G(z, L) can be shown to be finite every where for z £ [0, L/2] except at z = 0. It is straightforward 
to isolate the divergent part of G R (z, L) near z = 0, and we find that 



G R (z,L) = A eM 2z 2z) +G' R (z,L) 7 



(43) 



where G' R (z, L) is finite for all z £ [0, L/2]. Hence, the first term of Eq. (|42|) is finite but the second has a logarithmic 
divergence at z = when A > 0. Introducing a large P cut-off in the integral in Eq. I|37fl defining G R (z, L), as was 
done in [9|, eliminates the 1/z singularity in G(L, z) but has no physical basis, in fact it is eliminating a physical effect 
which should be there due to the presence of image charges; the divergent term of G R (z,L) should be kept in the 
exponential (which is just the Boltzmann factor due to the image charges). We can now carry out a legal expansion 
in g for the terms which give finite contributions, and obtain 



/3 ^ 1) = 32^ G « (L/2 ' i) + l^ 



L/2 q l 

dz G' R (z,L)—G' R (z,L) + +— 



L/2 d 



1 — exp 



dL ' ' 4tt 

We now write the two-loop contribution to the disjoining pressure as the sum of these three terms 



gAe 



-2z 



4z 



(44) 
(45) 



The first two terms can be immediately written down and a careful asymptotic analysis of the last term of Eq. I|44|) 
gives a finite resummed two-loop contribution 



(la) _ 


9 


d " 


32tt 


,(16) . 


9 


d ~ 


8tt 



pp. 



P ' d " 167T 



f Q L/2 dz (g r (z, L)(-^G' r (z, L) - ^G^(0, L)) + G' R (z, L)-^-G' R (0, L) 



^G'zl ().!)( ln(^)+ I 1- 



(46) 
(47) 

(48) 



where 7 = 0.577215.... is Euler's constant. We see from the third term a contribution proportional to glng signifying 
that the assumption that Vd can be naively expanded in g is incorrect. 



III. LARGE L ASYMPTOTICS 



In this section we will analyse the asymptotic behavior of the disjoining pressure of thick films. We start by recalling 



the one loop result. At large L the behavior of the integral determining V d in Eq. (|40|l is dominated by small P 
the contribution to the integration coming from near P = 1. The leading term for large L is Q 



i.e. 



&vf „ _expJ_2L) m 

We see that the thick film disjoining pressure at one-loop is attractive and exponentially screened with characteristic 
length Id. Interestingly, the value of A does not appear in the leading term of this expression which is therefore 
independent of the strength of the dielectric discontinuity between the film and exterior. For large L we find that 

Gfl(L /2,L)«2^-^, (50) 
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and thus 



pp£ a) 



g exp(-2£) 
8tt L 2 



(51) 



which is clearly negligible with respect to the one-loop contribution at large L. 
The behavior of G'(z, L) at large L is given to leading order by 



G^(z,L)«exp(-2L) 



1 1 



2L- z 



exp(2z) 



1 1 



2L + z 



exp(— 2z) 



+ H(z) 



where 



Hence 



/oo 
dP (r(P) - A)exp(-2Pz) 



dL 



G' R (z,L) w -2exp(-2L) 



1 1 



2L- z 



exp(2z) 



1 1 



2 L + z 



exp(— 2z) 



Using this result we find that for large L the leading behavior is 



fmM ~ 9 — ex P(~ 2L ) 

M d ~ 16tt L 



ln(i) + 7 + 



1 — w 



(ln(l + n) -nln(2)) + 



1 HL) \ 
L 1 / 2 ' L / 



and for the special case n = the non-leading Ofd/L 1 / 2 ) correction is absent and this result becomes 



(ib) 



16tt 



exp(-2L) 
L 



ln(L) + 7 - 



2L 



ln(L) + O 



(52) 



(53) 



(54) 



(55) 



(56) 



It can be shown that the coefficient of exp(— 2L) ln(L)/L is independent of n in the general case and is therefore given 

for any value of n by this result. Hence, as in the case of the leading contribution to at large L is also 

independent of A. 
We also find 



PP, 



(lc) 



gA exp(-2L) 



[-ln( 5 A/2) + l-2 7 )] 



4tt L 

We notice that if we take the contribution of the one-loop term with the leading two-loop term then we have 



(TP 



which suggests the speculative resummation 



exp(-2L) 



4ttL 



1 
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ln(L) 



PP 



exp(-2L) 
" 4tt£ 1 -3/ 4 



(57) 



(58) 



(59) 



If valid, the terms leading to this resummation must come from higher-order terms in the cumulant expansion. It can 
be verified that they do not arise from higher-order terms in g in the expansion of Eq. (|41Jl . 

Examining the contributions to (3V ty1 ^ 1 as defined by Eq. (j45(l at large L we see that (3V^ la ^ is clearly subdominant, 
but although (3V^ lb ^ is clearly the leading term at very large L we must bear in mind that this expansion is valid for 
small g. Comparing Eq. (56) with Eq. (57) we see that there exists a cross-over length L* defined by 



(60) 



such that for 1 <C L* <C L we have (3V^ ~ /3"p( lc ) and thus in this region pP^ has the same functional form a the 
one-loop result. Thus only for L^$> L* does one see the modification of the functional form of pP with respect to the 
one-loop result. 
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IV. NUMERICAL RESULTS 



The major result of this paper for the large- L asymptotic behaviour of Vd is the expression Eq. (|55j) . In Fig. 
we show data points for 

R W {L) = ^e^/SpW (61) 

versus L computed numerically using Eqs. (|37|) . 143|) and l|47|l for n = 0. The solid line is the asymptotic form derived 
for R d lb) {L) from the expression in Eq. I|56|) but including a term proportional to 1/L whose coefficient we determine 
by a fit to the data. This fit curve is 

Rf\L) = - 7 -h(I)-ih(I)-if*. (62) 

It is not useful to directly plot given in Eq. I|41l) since it is a rapidly changing function of L. However, the 
important comparison is with the one- loop result vP , which corresponds to the attractive Casimir force and which is 

present even as p — > 0, that is, g — > 0. In this limit vP arises solely from the thermal fluctuations of the electrostatic 
field in the presence of interfaces of dielectric discontinuity. In the limit g — ► the two-loop contribution must vanish 
since the one-loop result is exact; this can be seen in the expressions given above. However, for g > we see that 
l^P is of the same sign as vf" 1 and so corresponds to an increase in the attractive force. In Fig. we plot the ratio 

^rf 1 ' V^d f° r various values of g: g — 0.1,0.3,0.5. We see that V d is comparable with V d °^ for g ~ 0.5 over a wide 
range of values of the film thickness L. For small enough L (L < 2) the one-loop result eventually dominates. It is 
clear that VP scales, as expected, approximately linearly with g in the small range considered. The value g ~ 0.5 
corresponds to a solute density in water of about 50 mM at room temperature. This result shows that a quantitative 
analysis of the collapse phenomenon in a thin electrolyte soap film must take account of higher-order effects for solute 
densities greater than 10 mM. An important point to note is that both 7>P and VP are negative for all L > 0, and 
so the two-loop contribution adds to the attractive force between the interfaces. 

In Fig. © we plot the ratio V d /V d lc ^ where the denominator is the approximation to pfp given in Eq. (|57jl . This 
approximation arises from the resummed terms which properly account for the Boltzmann factor associated with the 
image charges. It is clear that this approximation is a good one over a wide range of film thickness L and shows that 
the dominant contribution in this range behaves as gln(g). It is clear that this should be the case in this range as in 
the case g = 0.5, from Eq. L* ~ 16. 



V. CONCLUSION 

In this paper we have carefully calculated the two- loop contribution, V d , to the disjoining pressure for a thin 
electrolyte soap film of thickness I and monovalent solute density p. The surfaces of the film are the interfaces of 
dielectric discontinuity separating the film interior of dielectric constant e from the exterior which consists of non-ionic 
medium of dielectric constant eo- The length scales in the problem are set by the Bjerrum length Is = e 2 (3 / (Aire) and 
the Debye length lp = m = y / 2pe 2 (3/e. The model used for the ionic interactions is the Sine-Gordon field theory 
extensively discussed in earlier work [8j . All quantities can be expressed in dimcnsionless form in terms of the coupling 
constant g — Ib/Id and n — eo/e, with overall dimensions carried by the Debye mass m. By naive power counting 
the N— loop contribution is the 0(g N ~ 1 ). The self-energy divergences are included in the explicit renormalization 
constant Z(g) defined in Eq. JSJ. 

We have computed V^p in the cumulant expansion of the film free energy ln(S^?) given in terms of the effective 
action Sf(4>,L) defined in Eq. (|12fl . Although, we are naively working to O(g) it is clear that it is illegal to expand 
all terms to this order since the effect of image charges is to introduce a spurious Ultra- Violet divergence in this case. 
As an ion at distance z from an interface approaches the interface the image potential is repulsive and diverges like 
(iVi(z) fa gA/2z, where A = (1 — n)/(l + n). It is necessary to keep the Boltzmann factor exp(— (iVi(z)) intact for 

this otherwise divergent contribution, which is then regulated in the proper way. The resulting expression for 

is given in Eq. (|44J) . The outcome is that the naive expectation that this contribution is O(g) is incorrect and the 

asymptotic analysis of the expression for V^p gives the term V d lc ^ shown in Eq. (|48|l which contains a term of order 

g\n(g). Moreover, in Fig. J2J we see that the approximate expression V d lc ^ contributes the major part of the full 
result for small L = I /Id, and so the collapse phenomenon observed in thin electrolytic soap films may be modified 



significantly by such non-analytic terms. From Fig. J5J) we observe that the two-loop and one-loop contributions are 
comparable for g ~ 0.5, corresponding to solute densities p ~ bOmM . For L< 1 we would still expect the one-loop 

contribution, to dominate all others on physical grounds; the image charges will ensure that all ions are expelled 
from the film leaving a non-ionic water-filled film. The only, and therefore exact, contribution is the Casimir term 
due purely to the presence of the dielectric discontinuities. 

The question arises how large g can be taken before the loop expansion breaks down. The only evidence is from 
the calculation of the bulk Debye pressure where the Debye-Hiickel correction to the ideal value Pq is 

P B = Po(l-f) • (63) 
6 

This suggests that the expansion parameter might be as small as <?/6. However, we know that the two-loop term is 
0(g\n(g)) and that our calculation deals with the effects of an interface whose physics is completely different from 
that of the bulk. Nevertheless, we may take this result as an indicative of what we might expect. 

We might have expected that the leading large- L contribution from V^P behaves in a similar manner to that of the 
one-loop contribution, namely, 

7>< 0) « -^exp(-2L). (64) 

However, the analysis of the large L behaviour shows a departure from expectations because of subtle contribution 
from the image charges; we see from Eq. I|56|) and from Fig (JTJ that the leading large-L term behaves as 

Vf ] « -gA ,] ^-eiq>(-2L) . (65) 

This result implies that for sufficiently large L the loop-expansion fails to converge since will become larger than 

. We suggest that higher-order loop contribution can be resummed to give an anomalous exponent as shown in 
Eq. 159JI : this point requires further investigation. However, for physically relevant conditions the loop expansion will 
be valid. 

The conclusion of this paper is that the naive loop expansion of the free energy for these systems is invalid in 
the presence of a dielectric discontinuity, and a careful resummation of higher loop-order terms must be employed 
to properly take into account the Boltzmann weight of image charges. Furthermore, at two-loop order as generated 
by the cumulant expansion, the disjoining pressure is finite without the need for the introduction of an Ultra- Violet 
regulator. In addition, we have shown that the two-loop disjoining pressure is negative for L > and so enhances the 
attraction between the interfaces predicted at one-loop order. 
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FIG. 2: The ratio of the two- loop contribution, and the one-loop contribution, Vf\ to the disjoining pressure is plotted 

for different values of g — Id/Ib- g = 0.1, 0.3, 0.5. We see that is comparable with Tf^ for g ~ 0.5 over a wide range of 
values of the film thickness L. For small enough L (L < 2) the one-loop result eventually dominates. It is clear that scales, 
as expected, approximately linearly with g in the small range considered. The value g ~ 0.5 corresponds to a solute density of 
about 50 mM at room temperature. 
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FIG. 3: The ratio of the exact two-loop result, for the disjoining pressure and the approximate result given in Eq. 

1571 . V^ c ^ is seen to be a good approximation until L becomes small or larger than L* as defined by Eq. 160H (L* ~ 16 for 
5 = 0.5). 
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